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STJMARY 


In this paper further results of a study of the neutral hydrogen 
distribution in the upper atmosphere are presented. These are based on a 
Monte Carlo calculation of the quasi -stationary neutral hydrogen dis- 
tribution in the transition region between the collision dominated 
lower thermosphere and the collisionless exosphere. In this region the 
rapid loss of fast atoms perturbs the velocity distribution of hydrogen, 
in effect, cooling the hydrogen relative to the ambient main atmosphere, 
thereby decreasing the relative loss rate determined by escape. Conse- 
quently, under stationary conditions, the hydrogen concentration above 
the production region is significantly enhanced, leading to a higher 
hydrogen content in the atmosphere above 120 Km. Results of new calcu- 
lations of the hydrogen concentration as a function of altitude are 
presented, as well as the hydrogen content between 120 Km. and various 
altitudes, for various temperatures of the thermopause, between 1000 
and 2500°K. The latter reveal that the "total" hydrogen content of the 
atmosphere is increased by a factor of ~ 1.1 (at 2500°K) to ~.2.0 
(at 1000°). This points to a greater swing in hydrogen content between 
temperature extrema than that deduced in the convential manner. Both a 
higher hydrogen content and a greater amplitude of variation with 
temperature appear to came closer to an explanation of the experimental 
evidence (mainly Lyman - a radiation) . But they also emphasize further 
the need for a calculation taking into account convection of hydrogen 
between regions of varying concentration and the effect of the diurnal 
variation of temperature, which is too fast to permit the establishment 
of stationary conditions throughout the whole atmosphere. 



INTRODUCTION 


Though considerable efforts have been spent over the last decade on 
the elucidation of the problem of hydrogen in the earth's vicinity (see 
e.g. References 1-5, and the numerous studies cited there) attempts at 
its solution have met only partial success, and many related observed 
phenomena still elude a satisfactory quantitative explanation. 

Ever since the pioneering rocket experiments on Lyman a radiation 
by Friedman and co-workers (References 4,5) and the interpretation of 
their results by Johnson and others (References 6, 7), it had become clear 

that a size^Ple amount of neutral hydrogen, in the order of 10^ protons 

-2 

cm , is attached to the terrestrial globe, moving with it through the 
interplanetary medium. This telluric hydrogen corona absorbs, scatters 
and re-radiates the incoming solar Lyman a radiation, producing the 
Lyman-a glow observed at high altitudes, both by day and at night. 
Moreover, its presence is operative in determining the structure of those 
components of the radiation belts which are sensitive to interaction with 
neutral hydrogen, such as low energy protons with high pitch angles at 
relatively great geocentric distances from the equatorial plane 
(References 8,9). 

Though theories of the exosphere (References 2, 10-15) throw con- 
siderable light on the variation of the relative hydrogen concentration 
above 2500 Km, Bates and Patterson (Ref. l6), and Kockarts and Nicolet 
(References 17-20) were the first to present a detailed picture of the 
neutral hydrogen distribution throughout the whole atmosphere, 
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They recognised. that Toss 


for various temperatures of the thermopause. 

of hydrogen incurred from thermal escape leads to a serious departure 

from diffusive equilibrium in the thermosphere, thereby affecting its 

vertical distribution at all altitudes. This then results in a strong 

12 -2 

variation of the atmospheric hydrogen content, from ~ 5X10 atoms cm 
in a column between 1200 and 2000 km, at a thermopause temperature of 
1000°K, to ~ 1 .3X10^ atoms cm ^ at 2000°K. 

This significant advance in the treatment of the problem enabled 
Thomas and Donahue (Reference 1, References 21-25) to carry through the 
solution of the difficult radiative transfer problem involved in Lyman-a 
radiation. Though the results of these calculations are indicative of 
the basically correct approach of the authors quoted in the last paragraph, 
some difficulties were left unresolved which mar somewhat the validity and 
applicability of the hydrogen atmosphere model proposed by them. 

It provided no satisfactory quantitative explanation of the observed 
intensity variation with solar zenith angle, nor could the required 
optical depths, which are proportional to the hydrogen content above 120 Km, 
be satisfactorily reconciled with the atmospheric conditions, mainly 
temperature, prevailing when the Lyman a radiation was observed. 

Though these difficulties are related to the time -dependent 
asymmetry of the hydrogen corona, and the associated lateral flow 
(References 2k, 25), we first address ourselves to a problem which appears 
of basic importance even in the time-independent treatment of a symmetric 
hydrogen atmosphere. As Nicolet points out (Reference 20), one cannot 
accept a Maxwellian distribution as representing real conditions at the base 
of the exosphere. Qne might also add that the concept of a fixed, 
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temperature independent critical level at ~ 500 Km, underlying the 
calculations in References 16-20 might have to he discarded, in view of 
the considerable variation in the number density of the main constituents 
at this altitude, which may change by about an order of magnitude, say, 
between thermopause temperatures of 800 and 1200°K (Reference 26). 

In a previous paper (Reference 27), a method was described to 
evaluate the hydrogen velocity distribution in the vicinity of the critical 
level, which takes account of the non- equilibrium situation, brought about 
by thermal evaporation. Preliminary results for thermopause temperatures 
of 1500 and 2500°K were also presented, indicating that the relative escape 
rate of hydrogen is, indeed, reduced quite significantly. As pointed out 
there, the reduction is greater than deduced in earlier attempts (References 
2, 28, 29) at solving the escape problem, since at no point in our calcu- 
lations any particular form of the velocity distribution was introduced 
in the calculation. 

The basic idea underlying our method was the division of the upper 
atmosphere into three regions, and can best be visualized with the aid of 
the schematic model in Fig. 1: lowermost is the bulk of the thermosphere, 

where solution of the Chapman- Cowling diffusion equation (Reference 30) 
appears adequate in deducing the vertical distribution of hydrogen. Next 
we evaluate the evolvement of the quasi-steady-state hydrogen distribution 
in an extended transition region by means of the Monte Carlo method. The 
uppermost region is the exosphere, where again conventional methods may 
be applied. In a sense the two upper regions overlap partially, since we 
place the critical level, the exobase, ' some distance below the top of 
the transition region, where the Monte Carlo calculation is carried out. 
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This procedure then serves to take some account of the occurence of collisions 
in the exosphere, rare though they are, and eliminates spurious effects 
resulting from the sharp artificial boundary at the top. This boundary 
bears the character of a semi permeable membrane reflecting slow atoms, 
but passing those above escape velocity. For further details of the method, 
and the underlying shochastic model our previous paper may be consulted. 

In this present note additional results for thermopause temperatures of 
1000 and 2000°K are presented, as well as the results of a new calculation 
of the steady-state hydrogen atmosphere above 120 km, based on the Monte 
Carlo results. In these calculations all three methods described in the 
last paragraph were applied, each in its appropriate range. 

The Monte Carlo results in the transition region for thermopause 

temperatures of 1000 and 2000°K fit in well with the earlier results . 

Those for 2000°K yielded values of the significant parameters (effusion 

velocity, response time of the hydrogen distribution to abrupt changes in 

temperature, concentration gradients etc.) intermediate between those of the 

bracketing temperatures of 1500° and 2500°K. The results for 1000°K also 

appear to be consistent with these, though in the latter case the statistical 

2 

accuracy of the effusion velocity, ~ 2X10 , is rather poor. 

Again, as in the earlier results, interpretation of the dispersion 
of the various velocity components as effective temperatures of the 
hydrogen, displays the " anisotropic cooling " effect, observed there. 

This interpretation appears to be quite useful in the evaluation of 
the hydrogen distribution in the exosphere. The relative concentration in 
the exosphere derived from the raw Monte Carlo data on angular and velocity 
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distribution at the exobase differs only by at most 10 % from the values 
obtained, for an equivalent exosphere computed in the conventional manner 
from a Maxwellian distribution of hydrogen at the base, whose temperature 
is lower than that of the ambient oxygen, assumed in diffusive equilibrium. 

It appears that the relative concentration in the exosphere is reduced at 
large distances by a factor of 2-k, depending on temperature. In the light 
of the cooling effect discussed, this reduction follows from the reduced 
hydrogen scale height caused by the lower temperature. 

The absolute concentration in the exosphere, however, is determined by 
the concentration at the critical level which has to be obtained by deducing 
the concentration profile of hydrogen throughout the thermosphere. This 
has been done by solving the diffusion problem in this region in a manner 
similar to that first proposed by Mange (Reference 31), and applied since 
in the calculations of References 16-20. We have used the Harris and 
Priester model atmospheres (Reference 26) as the medium through which 
hydrogen diffuses upward. 

Finally, the three separate portions of the concentration profile were 
combined and the concentration profile of hydrogen above 120 Km was obtained, 
for the four thermopause temperatures considered. Integration of the 
concentration curves between 120 Km, and higher altitudes yields then 
the hydrogen content between these two layers . 

Since Harris and Priester (Reference 26) base their models on constant 
boundary conditions at 120 Km, while Kockarts and Nicole t (References 17-20) 
use constant boundary conditions at 100 Km, we have carried out corresponding 
calculations based on the diffusion velocities derived from Jeans' formula. 
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This permits comparison between the profiles obtained when the effect 
of escape on the hydrogen velocity distribution is taken into account, and 
those deduced in the convential manner. 

The important results from this comparison are as following: in the 

thermosphere the concentration of neutral hydrogen is effectively closer 
to diffusive equilibrium than convential escape theory predicts, while the 
converse holds in the exosphere, a consequence of the cooling effect 
discussed. Overall, because of the reduced relative escape rate, and the 
ensiling prolonged sojourn of hydrogen atoms in the atmosphere, concentrations 
are, in general, higher than obtained hitherto, except at large distances, 
in the order of 4 earth radii. Consequently, for the same rate of supply by 
the source of atomic hydrogen from photo-dissociation of water vapor and 
methane (Reference 31 ), the atmosphere has a higher hydrogen content, than 
predicted hitherto. 

6 -3 

If we assume an average concentration of 10 cm hydrogen atoms at 

2 

120 KM, the respective values of the hydrogen content in a 1cm column 

above this level are ~ 1.37X10"^ cm ^ at 1000°K and 2 . 00X1 0 1 ^ cm ^ 

o 12 12 ”2 

at 2500 K, as compared to 6. 9X10 cm and 1.79X10 cm at the 

corresponding temperatures, when the conventional loss rates are applied. 

More details can be inspected from the curves in the following sections. 

These results suggest higher optical depths for Lyman cl , and predict 
a greater swing of the hydrogen content than assumed up to now. Though this 
appears to be in the right direction (see e.g. Reference 1, References 21, 22), 
it even more stresses the need for a comprehensive time dependent treatment 
of the problem. 
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FURTHER RESULTS OF THE MONTE CARLO CALCULATION 
ON THE HYDROGEN DISTRIBUTION IN THE TRANSITION REGION 

In this section the results of additional Monte-Carlo computations 
of the hydrogen distribution in the transition region are represented, 
and some general conclusions are drawn from these results. Since the pro- 
cedure applied was identical to that reported in a previous communication 
(Reference 27), and so are the definitions of some relevant quantities, 
this should be referred to for details. 

As was done there, the atomospheric parameters used in the compu- 
tation are based mainly on the work of Harris and Priester (Reference 26), 
and are summarized in Table 1. 

TABLE 1 

ATMOSPHERIC PARAMETERS USED IN THE 
COMPUTATION FOR THE 1000° AND 2000°K 
CASES 


Atmospheric 

Temperature 

(*K) 

Altitide 

(Km) 

Oxygen 

Concentration 
(cm ^) 

Mean 

Free Path 
(Kin) 

Scale 

Height 

(Km) 

1000 

320 

0 

3 . 5 X 10 

10.54 




6 


62.5 


580 

5.5x10 

670.00 


2000 

500 

O 

2 . 32 X 10 

£ 

19.5 

112.0 


95 0 

k. 20 X 10 

1080.0 



Again, as in the previous work, the limits of physical validity, 
i.e. the boundaries of the physical transition region between thermosphere and 
exosphere, are taken at the lower end where the slope of the M. C. curve approa 
ches that of an approximate diffusive equilibrium curve, while at the top 
the exospheric base (A=2.5H) is assumed to bound it. 
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The part of the density curve obtained from the Monte Carlo computation 
in the transition region is represented by the beaded portion of the solid 
curves in Figures 1 and 2, extending between 300 and 480 Km., respectively 
520 and 800 Km for the two additional cases treated. As in the previously 
reported cases, also included in these figures, there is an inflexion in 
these curves, which reflects the closer approximation to diffusive equili- 
brium in the lower thermosphere, compared to conventional calculations, 
whereas toward the exosphere the departure from equilibrium is more pronounced. 

Figure 4 represents the response or relaxation times of the layers 
near the critical level. These can be seen to vary inversely to temperature 
between ~ 1300 sec. to 1000°K to ~ 65 O sec. at 2500°K. They can best be 
interpreted as characteristic times for the establishment of stationary 
conditions in the velocity distribution of hydrogen in the transition region, 
while the response of the concentration will be determined to a considerable 
extent by the altered vertical and lateral diffusive flux into the region. 

This effectively means that the change in relative effusion velocities, 

shown in Figure will lag behind the change in temperature approximately 
-t/T 

proportional to e ' * while the change in the absolute loss rate will vary 
in a more complicated way. 

Wow, with respect to the relative reduction in effusion velocity, 
compared to values derived via Jeans formula, it has alreay been observed 
in our previous paper that this reduction is greater for the lower thermopause 
temperatures, in the range considered, than for higher ones. Some reflection 
may convince us that this result is not quite so surprising as one might 
be tempted to think at first. 
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It is quite clear that Jeans 1 formula furnishes upper limit values for 
thermal escape (see e.g. Reference 33). Consider now the hypothetical case 
of zero temperature. Evidently there is no thermal escape whatsoever, 
neither as derived from Jeans 1 formula for E->°°, nor from a hypothetical 
Monte Carlo calculation, if this were possible. Thus we may assume that 
at very low temperatures the two curves for the effusion velocity approach 
each other. 

On the other hand consider a relatively high temperate T. Let us 
consider the maximum geocentric distance R attainable by hydrogen atoms 
of mass m and mean energy 3Kr/2k where k is the Boltzmann constant. From 
concerva' .on of energy we obtain 

3hT/2-mg(R ) R = -mg(Rm)Rm (l) 

o o 

= -mg(R ) R 2 / Rm 
o o 

where R is, say, the geocentric distance of the vertical level, and 
o 

g the gravitational acceleration. Then 

Rm = R E/(E-1.5) (2) 

where E is the escape parameter of (R ) R /kT. Equation (2) means 

o o 

that once E reaches a value of 1.5, eventually all hydrogen will escape 
from the atmosphere. Thus escape will be very fast, and the two effusion 
velocity curves will again intersect. This is illustrated by the 
extrapolation of the Monte Carlo curve beyond the computed value for 
2500 k. In between there is a region of maximum departure of the real 
effusion velocity from the Jeans effusion velocity. This seems to be 
located towards the lower thermopause temperatures in the range considered 
here. 
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Moreover, we may also consider the cooling effect discussed above. 

In a first approximation let us assume that the fractional decrease in 
hydrogen temperature, T eff, at the critical level is proportional to the 
ambient thermopause temperature, T®, i.e. 


T eff = T “ (l ' a) 


(3) 


If we apply Jeans' formula, we obtain for the ratio, R , between the 

* y* 


effective effusion velocity, and the Jeans effusion velocity 




R v = (l-a) [l+E/(l-a)j e 


aE 

lTa 


/ (l+E) 


(4) 


If we take ctr-0.2, we obtain for hydrogen in the terrestrial atmosphere the 
values of R^ given in Table 2, as a function of temperature. 

TABLE 2 


-v- 


4-M k/J. J. U1U UiillOXUU V OJ 

T (°K) 

E 

R 


V 

1000 

7.05 

.265 

1500 

4.64 

.423 

2000 

3-26 

.552 

2500 

2.77 

.606 

3000 

2.31 

.6 58 

4000 

1*73 

.729 

4500 

1.54 

.814 


Again we observe that the relative effect is more pronounced at lower 
temperatures, in agreement with our computational results. 

Figure 6 represents the effective temperature of hydrogen at the critical 
level as a function of the thermopause temperature. This is defined here 
as the geometric mean between the critical effective temperature and the 
omnidirectional effective temperature, defined in Reference 27 . The rationale 
in applying this particular averaging method was the desire to take some 
account of the geometry of the actual problem, whereas the computations 
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were carried out in plane geometry. After these averages were taken 


it was found that the exospheric distributions decuced with these tempera- 
tures yield an almost perfect fit to those obtained from the raw data of 
the Monte Carlo computation. 

It is seen that in the range of our computations the effective temperature 
lies about 25$ below the ambient thermopause temperature. We believe 
this effect to be real, and consider it to be the combined effect of both 
the rapid escape of fast particles, as well as some loss of speed by 
fast hydrogen atoms in collision with oxygen, since collisions are not 
frequent enough to bring about detailed balance . If there were some heating 
of the ambient oxygen it would be too minute to produce any observeable 
effects. The extrapolation of the effective temperature curve beyond 
the computed values is of a tentative character. Beyond T o ~4500°K (E~1.5) 
escape presumably becomes so fast that the concept of temperature may no longer 
be meaningful. 

Though no Figures are shown here, representing the "microscopic" 
results, such as angular and velocity distributions for 1000 and 2000°K, 
these bear the same character as those of the earlier results for 1500 and 
2500° which are amply illustrated in the earlier paper (Reference 27). 
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THE RELATIVE CONCENTRATION OF HYDROGEN IN THE 

EXOSPHERE 

The relative concentration of an atmospheric constituent in the 
exosphere may be evaluated, once its velocity distribution at the critical 
level is known. The latter is defined by flpik andSinger as the level where 

HnojO-5 (5) 

Here H denotes the real height of effective atmospheric number density 

(i.e. that of oxygen in our case); no is the effective number density at the 

leval, and a an effective collision cross section. 

A more precise averaging method used by us does not materially change 

this criterion. If the angular probability density of the upward flux is 

given by 2xdx (with x the zenith angle cosine), the condition for a collision 

probability of jjs- in the upward omnidirectional flux is 

1 

2 J e" H ^'' oX xdx=.5 (6) 

o 

changing to a variable |=H/X. , one obtains 

e' § -|e _§ -§ 2 Ei(-|)4 (7) 

then Ei(-§) is the exponential integral 

00 

-Ei(-S) = J - e_^ dt 

S t 

Equation (7) has a solution §~0.4 in fair agreement with Equation (5)* 
Applying this criterion we obtained with the aid of the atmospheric data 
used altitudes of 480, 580, 800, and 9<X> Km for the altitude of the critical 
level at respective thermopause temperatures of 1000, 1500, 2000 and 2^00 K. 
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In the Monte Carlo computation we obtained the required velocity 
distributions for hydrogen in the form of two dimensional arrays c ., 

j- j 

essentially histograms of the population of thin layers (5-10 Km thick. 

Here the first subscript, i, refers to one of the 20 equal intervals in 

(-1, l), the range of direction angle cosines, while the second subscript, 

j, refers to one of the 50 intervals in speed between 0 and 10kT<»/m, 

the latter varying from ~ 40 Km/sec at 1000°K to ~ 6k . 3 Km/sec at 2500°K. 

To facilitate computation in the exosphere the "vectors" g. . for any 

J 

particular value of i were numerically approximated by a least square 
method by functions 

j. (v) = V 2e ' SiV (8) 


A goodness of fit test was applied to the approximations thus obtained. 

As observed in our previous paper the distributions obtained are quite 

symmetrical (though not isotropic) in x, the zenith angle cosine, up to voo. 

The relative concentration at an inverse geocentric distance y=R /R may 

o 

then be represented by a sum of integrals (similar to the double integral 
in References 11 or 14) 

P(r) = y 2 \ F i (n ± , y ) (9) 


where 


P(y ) = (R)/ n^ (R ) 


“i ■ <x i > 


a mean value of x in the ith interval, and 
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dv 


°° f r 1 

F iW y) = j fi K ) V 2 [i -f (i-m-J) ]-^(i-y)t ^ 


v min (n , y) 
2 


fi(v) {v 2 Ll-y 2 (l-q 2 ) ]-v 2 (l-y)} v 5 d 


( 10 ) 


v min ((i.jy) 

v (li^,v) is obtained from conservation of energy in orbital motion, as 
min H 2 


v min = V » (l_y)/[l ' y (1 ^ 2)1 


( 11 ) 


To test the dependence of tbe exospheric distribution on the exact 

choice of the escape level - and this is important in our study since we 

do not use an exact Maxwellian distribution there - we carried out the 

calculation of Equation ( 9 ) with the values for F^ obtained from the 

arrays g.. . in two nearby layers, ~ 50 Km apart, near the "exact" escape 
J 

level. The resultant exospheric concentrations were almost identical 
everywhere within ~ 3$ of each other. 

In a similar manner exospheric concentrations were derived by using 
at the escape level Maxwellian distributions appropriate to the effective 
temperatures represented in Figure 6, which also closely approximated 
the curves for the Monte Carlo derived exosphere. 

These curves are shown in Figures 7 and 8. Figure 7 illustrates 
the decrease in relative hydrogen concentration versus R/Re, where Re 
is the earth's radius. The dashed curves represent the variation in 
number density obtained from Maxwellian distributions at the respective 
ambient thermopause temperatures. To complete the picture, curves for 
the corresponding barometric distributions are inserted. As can be seen, 
the difference between the set of solid curves and the dashed ones increases 
quite rapidly with increasing geocentric distance; but the ratio between 
the concentrations approaches fairly constant value s at around 10 earth radii . 
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This can be seen more clearly in Figure 8, where this ratio R q is plotted 
versus S = R-R^/H^, which appears a more suitable parameter for comparison 
of results for different temperatures. Here H^ is the hydrogen scale 
height . 

As with the effusion velocity, the effect on the exosphere is more 

pronounced at lower temperature, the ratio R^ at large distances varying 

from about 0.25 for 1000°K to about 0.6 at 2500°K. Again the explanation 

is that a roughly constant fractional decrease in the temperature at the 

_ ( R-R ) /I! 

critical level is more effective due to the exponential factor e ' o', 

which enters into most calculations. 

We nave included Figure ( 9 ) showing the passage from the transition 
region to the exosphere. These curves represent relative hydrogen 
concentration versus S, normalized to unity at the assumed escape level. 

The beaded portion of the curve is obtained from the Monte Carlo calculation, 
while the solid portion shows the exospheric results discussed in this 
section. The region of overlap has been mentioned before, and has no physical 
significance beyond partial consideration of the rare encounters in the 
exosphere . 
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THE ABSOLUTE CONCENTRATION OF HYDROGEN ABOVE 120 Km. 


A knowledge of the relative variation of concentration in the transi- 
tion region and the exosphere is of limited practical value unless the 
absolute concentration at some level can be estimated. To do this, 
one has to calculate the concentration profile through out the thermo- 
sphere, down to the vicinity of the production region, where a more 
or less reliable estimate of absolute concentration can be attempted 
(References 17, 32). 

In doing so we proceeded in a manner similar to that applied by 
Bates and Patterson (Reference 16), starting with the equation of 
continuity 


n(R)w(R)R 2 =S R 2 


a a 


( 12 ) 


Here w denotes the net flow velocity of hydrogen at R S is the specific 
loss rate of hydrogen at a level , considered in the following as the 
upper boundary for our computation. We choose R at a level within the 
transition region, where the slope of the concentration profile approaches 
diffusive equilibrium. From there downward we feel justified in using 
Cha pman ' s expression for the diffusion velocity (Reference 30, Equation 14.1.1) 


w = - 


Vo 


^ j"d , n h x j_ V^j ( \ d log P j v 1 3T | 

5 115 + — ( vV “SR + -Tf 35 ; (13) ~ 

ho 


Here n^cm J is the hydrogen concentration, n o cm stands for the concen- 
tration of main atmospheric constituents - mainly atomic oxygen close to 
the transition region, and molecular nitrogen near 120 Km. 

Dgo is the binary diffusion coefficient for the 0-H mixture, 

nqu + n , m and m are the respective masses of hydrogen and an average over 
no o 

the main constituents, p is the pressure, T the temperature, is the 
thermal diffusion ratio. 
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In this scalar form Equation 13 accounts only for radial flow, 
corresponding to a spherically symmetric steady state model of the 
atmosphere. Now since n^^n^. Equation 13 reduces to 

% 


w = -n 


Ho 


1 ^4- <Hog p m 1 dp , „ dT 1 
\ dE - n ffi “dE “ r p dE + dE J ilk) 


where a T is the thermal diffusion factor, a T ~ K T n/n^ 
Since p = nKT 


dlog p 1 dn + 1 dT 

dE n dE T dE 


whence 


t.t - _r» -fi_ \ d log T m dlog p | 

W D ED \ dE U dE ~ J 


dE 


(15) 


(16) 


(17) 


Let z = Ep-E (E<Ep) such that Equation (12) becomes 

2a ^a 2 

w(B) = n(B) (Ea - z) 

This may now be combined, with Equation (17) to yield 


m 


E 


[§- + (l+a T ) d T " m K~ 1j % ^ j 

dz dz n D 12 U) v Ea-z 

where H is the local pressure scale height 

Z 

which may be used for an integration factor to obtain a solution of 
Equation (19) in the form 

z 

n^z) =n^(0){(T o /T z ) (l+aT) [1 + | 0(z 1 , e P(z ’ } dz 1 ]} 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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where 


and 




(23) 

(24) 


D„ was approximated by the hard sphere formula 
Mo 


D S = I n o U) " ° 2 £ 


tt 


numerical evaluation of Equation (21 ) yields the concentration of hydrogen 

below R , relative to that at R .. In the actual numerical work use was 
a a 

made of Harris and Priester's model atmosphere (Reference 26). Specifically 


the following tabulations were applied at the respective thennopause 


temperatures : 

lOOO^K Model 5, 00 hrs. 

1500°K Model 8, 08 hrs. 

2000°K Model 10, 12 hrs. 

2500°K extrapolation 

If one considers the relative concentration, assuming n ( R a ) =1 , S & becomes 

numerically equal to the net flow velocity at R & . To permit comparison 

Equation 21 was integrated twice for each thermopause temperature, once 

using a value S ^(M.O)based on the Monte Carlo result, while for the second 
0# 

calculation a value S (M.B.) derived from Jeans* formula was substituted. 

a 

Table 3 presents the respective values of Z , S , R and R , the assumed 

O Q» o 

altitudes of the critical levels for the various temperatures 
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TABLE 5 



VALUES OF R 
a 

, S , AND R 
a' 0 

FOR VARIOUS TEMPERATURES T 

00 

Too(°K) 

R (Km) 
a 

Z (Kin) 
0 

S (M.c.)(cm 2 sec 
a 

_1 ) S (M.B.)cm‘ 2 sec 
a 

1000 

6770 

k 8 o 

i.8oxio 2 

7.56X10 2 

1500 

6830 

580 

1. 85x10^ 

G . k 5X10T 

2000 

6970 

800 

5.05x10^ 

1.70x10? 

2500 

6970 

900 

1.20x10 

2 . 99X10 

For a T a 

mean value of 

-0.2 was taken from Reference 33. 

Equation 21 


was calculated from the respective level R downward to 120 Km, where 

a J 

the concentration of the various atmospheric constituents is constant. 


in accordance with the Harris-Eriester model (Reference 16). 

Since we are interested in absolute concentrations, the hydrogen 
concentration at this level was assumed to be 10^ cm Estimating the 

hydrogen diffusion velocity there as ~ 2.5 cei see ^ this gives a source 

7 -2 

flux of ~ 2.5X10 cm sec x , in agreement with Nicolet's earlier estimates 
(Refs. 17, 31, and 32). The three independent portions of the concentration 
profile were then combined and are presented in Figures 2, 3, 10 and 11. 

Solid curves are based on the M.C. results, the portion with solid circles 
being obtained directly from the Monte Carlo calculation in the transition 
layer, whereas the portion with open circles represents the exospheric 
distribution based on the former. Dashed curves represent the concentration 
profile based on the conventional treatment of escape. 

As a comparison of the two sets of curves indicates, what we believe 
to be the real hydrogen concentration profile above 120 Km is actually 
overall closer to diffusive equilibrium than a profile derived with the aid 
of Jeans' formula . This strengthens our contention that on a large scale escape 
perturbs the medium less than application of the former indicates. 
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Consistently the Monte Carlo concentration curves lie higher than 
the corresponding alternate curves, except at large distances (~ 4Rp) 
where the curves approach each other or even intersect. This is easier to 
accept than the converse; loss from escape should affect the concentration 
more at high altitudes than in the thermosphere at large kinetic depths 
telow the escape level. 

Inspection of the curves reveals that the difference between the 
corresponding concentrations reaches a maximum in the vicinity of the 
critical level, where the assumed sink is located. The ratio of the 
corresponding concentrations is in the order of closely related to 
the ratio of the corresponding effusion velocities. 

Of greater interest in direct relation to some observed phenomena 
than the concentration itself is its integral, J n(z)dz, between two 


specific limits, yielding the total hydrogen content between these two 
layers. Numerical integration of n(z) has been performed by us, the 
resulting curves being shown in Figures and 


N 


These confirm our expectation that the hydrogen content per cm , 
z 

(z) r 

' ' n(z' )dz', between 120 Km and higher levels (z) is considerably 

^ J 

e nha ncer! due to the smaller relative loss rate, as compared to results 
obtained in the conventional manner. 

Both in absolute numbers, and relatively, this enhancement is more 

pronounced at lower thermopause temperatures than at higher ones. This is 

00 

indicated in Table 4, where the "total hydrogen content per cm^, fn(z)dz 

120 

is presented both from the Monte Carlo results .as well as from the 


conventional calculation. Moreover their ratio R is given, and the ratio 
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Noo(T)/Nco(lOOO) of the total content at a thermopause temperature T°K 
to that at 1000°K. 

TABLE 4 

TOTAL HYDROGEN CONTENT cm -2 AND VARIOUS RATIOS FOR DIFFERENT 

THERMOPAUSE TEMPERATURES 



Nco(T) 

. -2 
i cm 


Nca(T)/Noo(lOOO 

T®( 0 K) 

M.C. 

M.B. 

R 

n 

M.C. 

M.B. 

1000° 

1.34X10 15 

6 . 9X1 0 15 

1.94- 

1.00 

1.00 

o 

1500 

4. 13X1 o 12 

2.36X10 12 

1-75 

0.309 

0.342 

2000° 

2.60X10 12 

1.89X10 12 

1.38 

0.194 

0.274 

2500 ° 

J.2 

2.00X10 

1.78X10 12 

1.12 

0.149 

0.258 


Of special interest in these results is the greater amplitude in 

_p 

the variation of the total hydrogen content cm with temperature. 

Though it cannot he assumed that the full swing of this variation can he 
realized on a diurnal scale, as will he discussed in the last section 
of our paper, it certainly predicts a larger difference in the hydrogen 
content between periods of high and low solar activity. This actually 
may he even greater than the above numbers indicate, since the latter 
are based on constant boundary conditions at 120 Km, while the actual 
situation may allow for some variation of hydrogen concentration even 
at this low level (see e.g. Reference l8). 
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CONCLUSIONS 


The results of this new investigation of the steady state neutral 
hydrogen distribution above 120 Km presented here indicate that thermal 
escape is less effective in disturbing the diffusive equilibrium of 
hydrogen in the thermosphere than was previously assumed, whereas, on 
the other hand, its effect on the distribution of hydrogen in the exosphere 
is stronger than thought hitherto. 

Consequently, in view of the smaller relative escape rate the overall 
hydrogen content of the atmosphere is enhanced by a factor of ~ 1. 1-2.0 
in the range of thermopause temperatures considered here. From the trend 
indicated in the variation of this factor it is likely that the relative 
enhancement may be in the order of ~ 2.5 for the range of temperature 
prevailing near solar minimum, which appears to be about 6oo-800°K. 

It should be noted, however, that the difference between the present 
model and previous ones is not great enough to permit on the basis of 
available experimental data an unequivocal decision with respect to its 
fundamental correctness. 

Indications of evidence in its favor appear to be the following: The 

present results can more easily provide the required optical depths called 

for a satisfactory explanation of Lyman ~ a experiments. (References 21-25, 

15 

References 54, 55)* These require the presence of up to ~ 4X10 hydrogen 

-2 

atoms cm above 120 Km, which can be well accounted for by our model, 

7 -2 -1 

without stretching the source requirements unduly beyond 2.5X10 cm sec 
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It appears easier, for instance, to deduce from our model the con- 
centrations used by Donahue (Reference 23) in his interpretation of the 
Javelin rocket experiment (Reference 36), carried out in January i960. 

li -]± ll 

These were ^ 9X10 cm and 5X10 cm J at altitudes of 500 and 1100 Km 

respectively. The thermopause temperature assumed by Donahue was 1000°K, , 

which appears 10-15$ too low (e.g. Ref. 37) > such that the concentration 

7 

at 1000 Km would have to be increased even above the high value of 3X10 

-3 

cm quoted by him. To yield the required concentrations, at the higher 

7 -3 

altitudes, 2X10 cm hydrogen atoms at the 100 Km level would suffice 
in our model. 

Nonetheless, our steady state model suffers from the inherent defect 

of all such models. In view of the even wider swing in the total hydrogen 

content with varying temperature it is hard to reconcile the required supply 

of hydrogen in sufficiently short time ( 10 hours ) to make up for the 

vast difference in the steady state abundance between day time and night 

time temperatures. This is so, notwithstanding the fact that the relative 

distribution in the transition region displays relatively very short response 

times, in the order of an hour or less. But the difference between the 

content cm” 2 between, say 1500° and 1100° (corresponding to conditions 

12 -2 

during the winter of 1959"60) is in the order of 5X10 cm , while the 

7 -2 -1 

assumed source yield is ~ 2.5X10' cm sec half a day. This is short of 
the required quantily by almost an order of magnitude. Consequently the 
full swing between the steady state abundances corresponding to diurnal 
temperature extrema will never be fully realized. 
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This has been recognized by Hanson et al (Reference 24), in their 
study of diurnal variations (Reference 24). Their conclusion that the 
day time temperature determines the mean diurnal concentration appears, 
however, questionable in view of the observed wide variation in day time 
concentrations (Reference 35 )> 'the relative variation of hydrogen abun- 
dance is greater at lower temperatures. 

Though their treatment of lateral flow was an important start, con- 
clusions from their and later attempts by Donahue and McAfee (Reference 25) 
have not yet sufficiently clarified to what extent the lateral diffusion 
of hydrogen is effective in limiting the diurnal variation. It appears that 
no satisfactory picture of the hydrogen profile in the geocorona can be obtained, 
before a time dependent, three dimensional model is investigated without 
neglect of transport in the thermosphere. 

Such an investigation may, for instance, throw light on the following 

effect: While transport in the thermosphere and near the critical level 

-2 

tends to reduce the diurnal difference in hydrogen content cm , this 
reduction may be limited in the exosphere. At some distance from the 
earth the nocturnal exosphere derives actually from day time hydrogen 
near the critical level (and vice versa). As a result of the latter's 
higher temperature it falls of less sharply than would be expected from 
night time temperatures . 
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Figure 1. Schematic Modal of tho Upper Atmosphere 
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Figure 2. n Hydregen Concentration (cm - ^) between 120 and 1500 km 
far 1000°K and V2000°K. Solid Curran- Baaad on M.C. Results. 
Dashed Curras- Conventional Calculatien. 
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Figure 3. a H , Hydregen Concantratiaa between 120 and 1^00 kn for 
T^» l£fc>0°K and T„- 2^00®K. Solid Currea- Baaed oa M.C. Reeulta. 
Daahed Curvaa- Conveatienal Calculatiea. 
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Figure 7. n/n Q , Relative Hydrogen Concentration in the Exosphere vs. 
Geocentric Di stare* ( in Units of Earth Radii) for Various Values of 
T . Numbers 1,2, 3, U Refer t* 1000®, 1500°, 2000°, 2500°K Respectively. 
Letters a,b,c Refer t* the Barometric Formula, the Conventional Exo- 
spher and an Exe sphere based on the M.C. Results, 




Figure 8. R n , Ratio of Exospheric Hydrogen Concentration* Derived from 
M.C. Results to Concentrations Derived from Conventional Theetf, ▼». 
s*=(r-r 0 )/H 0 , for Various Too, 







Figuro ID. 1m, Hydrogoa Concentration (cm“^) vs Altitude above 120 km 
for TJ1000®K and !1£2000®K. Solid Curvea- Baaed on M.C. Reaulta. Dashed 
Curvea- ConTontional Calculation* 
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Figure 11. riu, Hydrogen Concentration (cm“^) vs. Altitudo above 120 km, 
for T„» 1500°K and T,* 2500°K. Solid Curvoo , Baood on M.C. 
Results. Dashed Curves- Conventional Calculation. 
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